Yoctosecond photon pulses from quark-gluon plasmas 
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Present ultra-fast laser optics is at the frontier between atto- and zeptosecond photon pulses, 
giving rise to unprecedented applications. We show that high-energetic photon pulses down to 
the yoctosecond timescale can be produced in heavy ion collisions. We focus on photons produced 
during the initial phase of the expanding quark-gluon plasma. We study how the time evolution and 
properties of the plasma may influence the duration and shape of the photon pulse. Prospects for 
achieving double peak structures suitable for pump-probe experiments at the yoctosecond timescale 
are discussed. 

PACS numbers: 25.75.Cj, 42.65.Re, 12.38.Mh, 78.47.J- 



Recent advances in ultrafast laser optics aim at the 
creation of ultra-short pulses beyond the visible spectral 
range. Shorter pulses are desirable, since they provide 
better temporal resolution, while higher photon energy 
gives rise to better spatial resolution. High-order har- 
monics of femtosecond laser radiation have been shown 
to be sources of trains of attosecond extreme-ultraviolet 
pulses [J, Q] that can be used to produce single attosec- 
ond soft X-ray pulses [3, 3- For example, such single 
attosecond X-ray bursts [5j] have applications in molec- 
ular imaging 0,01; quantum control @, @|, or Raman 
spectroscopy (lOj. By introducing a controlled delay be- 
tween two such peaks, the dynamics of electron systems 
could be studied using pump-probe techniques Q. Such 
techniques also allow for the direct time resolution of 
many-body dynamics, like the observation of the dressing 
process of charged particles [ll| . In view of the obvious 
desire for even shorter pulses with higher photon energies, 
it is reasonable to search for alternative production meth- 
ods. In this spirit, it has been suggested that zeptosecond 
pulses could be created by focusing intense laser pulses 
on subwavelength-size structures [12j , by the reflection of 
a relativistically intense femtosecond laser pulse from the 
oscillating boundary of an overdense plasma or via 
nonlinear Thomson backscattering 14 1. 



Among the shortest possible time scales that are avail- 
able experimentally are those obtained through high- 
energy collisions. Particularly interesting in this context 
are heavy ion collisions that can produce a quark-gluon 
plasma (QGP), because they expose a complex system 
with rich internal dynamics that lives on a very short 
timescale. Heavy-ion collisions at the Relativistic Heavy 
Ion Collider (RHIC) and soon at the CERN LHC pro- 
duce this new state of matter up to the size of a nu- 
cleus (~ 15 fm) for a duration of a few tens of yoctosec- 
onds (lys = 10~ 24 s r* 0.3fm/c). In such a collision, 
the plasma is produced initially in a very anisotropic 
state, and reaches a hydrodynamic evolution through in- 
ternal interactions only after some thermalization time 
Tthcrm and isotropization time Ti so . The observed par- 
ticle spectra turned out to agree well with ideal hydro- 



dynamical model predictions [lR [ltl . which led to 
the assumption that isotropization times are as low as 
Tiso ~ 1 ys(0.3fm/c). However, it has been pointed out 
recently that viscous hydrodynamic models are still con- 
sistent with RHIC data if isotropization times as large as 
Tiso ~ 7ys(2fm/c) are assumed [3], even if the expan- 
sion before isotropization is assumed to be collisionless 
("free streaming"). Besiders a plethora of particles that 
are created in such collisions, also high-energetic photons 
are produced and detected experimentally figL | 



21|. 



So far, only time-integrated quantities have been consid- 
ered, like the extraction of the spatiotemporal extent of 
the emitting source through photon pair correlations [22 1 
and the suggestion to measure the isotropization time Ti so 
using the direct photon yield 23|, 24 1. 

In this Letter, we study the time-resolved production of 
high-energetic ultra-short photon pulses in the QGP. We 
demonstrate that the emission envelope depends strongly 
on the internal dynamics of the QGP. Under certain con- 
ditions, a double peak structure in the emission enve- 
lope can be observed, which could be the first source for 
pump-probe experiments at the yoctosecond timescale. 
We find that the delay between the peaks is directly re- 
lated to the isotropization time, and the relative height 
between the peaks can be shaped by varying photon en- 
ergy and emission angle. Such pulses could be utilized, 
for example, to resolve dynamics on the nuclear timescale 
such as that of baryon resonances As an alternative 
interpretation of our results, a time- resolved study of the 
emitted photons could provide a window to the internal 
QGP dynamics throughout its expansion. 

The QGP is formed in a collision of two heavy ions as 
illustrated in Fig. QJa-c) . We study the emission of direct 
photons from the expanding QGP 2fJ 2^1 . The energy 
spectrum of such photons extends to the GeV range, and 
the upper limit for the temporal duration of the GeV 
photon pulse is given by the expansion dynamics of the 
QGP, which leads to yoctosecond pulses. Regarding the 
shape of the photon pulse, we find that the angle-resolved 
photon emission rate strongly depends on the internal 
state of the plasma, characterized by the momentum dis- 
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Figure 1: (color online) Early stages of a high-energy collision, 
involving pre-equilibrium (first two columns) and equilibrated 
QGP phases (last column). Parts (a)-(c) show three snap- 
shots in time in position space. Shown are the two relativis- 
tically contracted colliding ions that create the quark-gluon 
plasma in the overlap region. Curly arrows denote photon 
emission and semicircles the detectors. Parts (d)-(f) are cor- 
responding pictorial representations of the plasma in momen- 
tum space. In an intermediate stage of the pre-equilibrium 
phase, the momentum distribution is anisotropic, resulting in 
a change in the angular photon emission pattern that can give 
rise to double-peaked photon pulses. 



tribution of the plasma constituents, see Fig. QJd)-[lJf). 
For an intermediate time after the collision, a momen- 
tum anisotropy occurs due to the longitudinal expansion 
of the plasma, which leads to a preferential emission of 
photons perpendicular to the beam axis z. A photon 
detector placed towards the beam axis would therefore 
measure a time-dependent photon flux. Based on a re- 
cent model il, 2fl for the internal plasma dynamics, we 
find that this mechanism can give rise to double-peaked 
pulse envelopes that could be utilized for pump-probe 
experiments. 

We base our model on the one-dimensional expansion 
model by Bjorken This assumes boost-invariant evo- 
lution of the quark-gluon plasma in a central region of 
the collision. The detector is placed away from the beam 
axis by an angle 9 within the reaction plane. Since we 
want to study hard photons produced in the GeV en- 
ergy range, we restrict ourselves to the production of di- 
rect photons emitted from the pre-equilibrium and equi- 
librated phases of the QGP. We will not consider photons 
produced in later stages of the collision where on average 
lower energy photons are emitted. For the same rea- 
son, we will omit higher order soft scattering processes 
like bremsstrahlung or inelastic pair annihilation, whose 
contributions are again more important at lower ener- 
gies 
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l23l |. The leading contribution to the photon 
production rate R originates from quark-gluon Comp- 
ton scattering and quark-antiquark annihilation. For 
anisotropic momentum distributions, the photon produc- 
tion rate R has to be calculated numerically [21(. This 



temperature T of the medium, the photon energy E 
and momentum k, the fine structure constant a, and 
the corresponding quantity for the strong force a s (we 
use h = c = k B = 1). The rate further depends 
on the anisotropy, which is described by a parameter 
£ = (Pt) I (pi)) ~ 1 that relates the mean longitu- 
dinal and transverse momenta ph and px (HI, I28I l29l|. 
To integrate this rate over time, we use a recent model 
to describe the time evolution of the pre-equilibrium and 
equilibrated QGP 3, H| ■ This model specifies the time 
evolution for the energy density £ = £(t), for the hard 
scale Phard = Piiard^) (which corresponds to T in the 
isotropic case) , and for the anisotropy parameter £ = £(r) 
as a function of the proper time r. These quantities scale 
with different powers depending on the QGP expansion 
dynamics. For example, the time evolution of £ can be 
written as £(r) = (t/tq) s — 1. In the free streaming phase, 
5 = 2, while in the ideal hydrodynamic phase 6 = 0. The 
model in Ref. [Sll] essentially introduces a smeared step- 
function for the exponent 5 = S(t) to interpolate between 
S = 2 at early times r <C T iso and S = at late times 
t 3> Tj so , where the duration of the transition is con- 
trolled by Ti so /7 with dimensionless parameter 7. In this 
way, both the pre-equilibrium phase and the equilibrated 
QGP phase of the expansion can be described by a single 
model. In the model, thermalization and isotropization 
happen concurrently, T th(!rm = r iso . 

In the radial direction, we use a straightforward gen- 
eralization of the temperature dependence T(r, tq) = 



T [2(1 - r 2 /i?r)] 1/4 [32, [33| for non-central collisions 



that is valid in the overlap region of the colliding nu- 
clei. As in Ref. [2^], we neglect the expansion of the 
QGP into transverse directions, since it is small in the 
initial stage throughout which the high-energetic pho- 
tons are emitted. To is the initial temperature, r is 
the distance from the center in radial direction, and 
Rt is the transverse radius of the nucleus. The num- 
ber of photons of frequency u>d that arrive at the detec- 
tor at time td in the laboratory system can be obtained 
by integrating the photon rate along all possible light 
paths. The world line of the photons is parametrized as 
z(t;t d ,x k ,y k ) = x k ki + y k k^ + k(t - t d + d), where k is 
the spatial direction of the light wave vector k = (u>, k), 
&i and &2 are two direction vectors orthogonal to k and 
to each other, and d is the distance to the detector. The 
signal fd that arrives at the detector at time td is given 
by integrating the photon emission rate along the photon 
world line f d (k,t d ,x k ,yk ) = J^ d dtf{k,t,z(t;t d 
where /(fc,i, z) is the Lorentz transformed photon emis- 
sion rate E d 3 R/d 3 k inside the QGP, and vanishing out- 
side. One finally obtains the photon rate per unit time 
and frequency interval with A d the detector area, 



d 2 N(u d ,td) 
dt d duj d 



8nd 2 



dx k dy k u>dfd(ud,td,x k ,yk) ■ (1) 



rate, obtained in the form Ed 3 R/d 3 k, depends on the For the concrete calculation we followed the param- 
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Figure 2: (color online) Photon emission rate as a function of 
detector time r. Solid blue lines show a large isotropization 
time Tiso = 6.7ys = 2fm/c with 8 = 2 (free streaming model) 
while dashed red lines correspond to a short isotropization 
time ri so = to = 0.3 ys = 0.088 fm/c. Parts (a) and (b) dis- 
play emission at midrapidity (9 = it/ 2) for a central collision 
with impact parameter 6 = 0. Parts (c)-(f) show double- 
peaked photon pulses obtained for b = 9.2 fm or 12.2 fm, and 
the vertical dotted line indicates the position of the larger 
Tiso = 6.7 ys. In parts (c) and (d), the detector direction is 
6 = tt/4, in (e) and (f), it is = tt/8. 

eters as given in We show results for LHC ener- 

gies, where the initial conditions are given by a formation 
time To = 0.3ys(= 0.088 fm/c), the initial temperature 
To = 845 MeV, and transverse radius of the lead ion 
Rt = 7.1 fm. The critical temperature, where the QGP 
ceases to exist, is taken as Tq = 160 MeV. We study the 
free-streaming model (S = 2) with parameter 7 = 2. As 
in Ref. [Sll] , the isotropization time is varied in the range 
Tiso = To to t; so = 6.7ys(= 2fm/c). Both possibilities 
are not yet ruled out by RHIC data. We use the model 
with fixed final multiplicity, where the initial conditions 
are adjusted as a function of Ti SO in order to result in the 
same final entropy as for Ti SO = To- 

Figures [2]( a) and[2jb) show the typical time evolution 
of the photon emission rate for a central collision with 
emission angle orthogonal to the beam axis (6 = 7r/2). 
We present results at 2 GeV (3 GeV) energy, where the 
photon production from the QGP at midrapidity is 3 to 4 
times as large as the production from the initial collisions, 
and roughly 6 times (50 times) as large as the production 
from the hadron gas j22|. The origin of the abscissa is 
the time when a photon emitted from the center of the 
collision arrives at the detector. Photons arriving earlier 
originate from a part of the QGP that is closer to the 
detector. The pulse shape is mainly determined by the 



geometry of the nucleus with radius 7.1 fm. 

Internal QGP dynamics occurs on a timescale of a few 
fm/c. Any structure of this order is washed out simply 
by the time for light to cross the size of the QGP. We 
apply two strategies to overcome this limit. First, we re- 
duce the physical extent of the QGP by considering non- 
central collisions with impact parameter b. Second, an 
optimization of the detection angle minimizes the trav- 
eling time through the plasma. In forward direction, the 
initial shape of the QGP is Lorentz-contracted, hence 
light travels through the initial shape quickly. This is 
partially spoiled due to the QGP expansion in the same 
direction. Thus intermediate emission angles are most 
promising for which the QGP appears partly Lorentz con- 
tracted but does not expand towards the detector. 

Figures E](c)-[2^f ) show the photon emission in a di- 
rection away from midrapidity 9 ^ ir/2 for a non-central 
collision. The impact parameter b = 9.2 fm (b = 12.2 fm) 
corresponds to an overlap region of the size 5 fm x 10.8 fm 
(2fm x 7.3 fm). Here, a striking new double-peak struc- 
ture appears. Roughly speaking, the minimum between 
the peaks corresponds to the maximum of the anisotropy 
parameter £(t), which grows with exponent 5 = 2 for 
t < Tiso, and vanishes for t > t; so . This is because the 
photon emission rate is suppressed for larger values of £ 
and smaller values of 9 [23j. The distance between the 
two peaks is therefore approximately governed by Tj so , 
indicated by the dotted line in Figs. [2lc)-[2jf). 

The dynamics shown in Fig. [2] requires a more de- 
tailed interpretation. In Figs. E]Jc)-[2jf) , the first peak 
corresponds to photons emitted from the blue-shifted ap- 
proaching part of the QGP, while the second peak corre- 
sponds to a slightly red-shifted and time-dilated receding 
tail of the plasma. The relative size of the first and sec- 
ond peak is controlled both by the emission angle and the 
photon energy. At increasing emission angles, the sec- 
ond peak is suppressed, because the phase space for pho- 
tons emitted from the receding part of the plasma is con- 
stricted. For smaller impact parameters or larger QGP 
sizes, the second peak overlaps with the first, thereby 
washing out this structure. Also for a short isotropiza- 
tion time Tiso = to (dashed lines) the separation into two 
peaks does not occur. Therefore this effect is very sensi- 
tive to the internal dynamics of the QGP. 

Our results are based on a number of model assump- 
tions. For example, we omitted the transverse expansion, 
which would become important at later times. It could 
either enhance the double peak structure through a fur- 
ther separation of red- and blue-shifted parts of the QGP, 
or reduce the effect through a prolonged photon pas- 
sage through the QGP. These modifications may cause 
qualitative changes to the photon emission envelope, but 
not to the yoctosecond timescale. In an actual exper- 
iment there is a background of photons from different 
sources [12]. These include photons produced by a jet 
passing through the QGP [33j, and could dominate the 
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effect that is expected from the QGP alone. Since these 
photons are produced on a similar yoctosecond timescale, 
they can be expected to have the positive effect of en- 
hancing the photon rate of the pulse. But at the same 
time they would constitute a background for the pulse 
shape determination. Other background photons are pro- 
duced at different time scales, e.g., due to decay of pions 
produced in the hadronization of the QGP, therefore they 
would not modify a time structure on the yoctosecond 
timescale. Photons produced from the initial collisions 
(pQCD photons) can be of comparable size, but they 
would only enhance the first peak of the double peaks 
depicted in Fig. E](c)-[2jf ) - A better theoretical under- 
standing of the QGP dynamics would be necessary for a 
more quantitative prediction of the expected pulse struc- 
ture. 

We now turn to an estimate of the photon production 
rate. In the GeV energy range, these are of the order 
of a few photons per collision [321]. Note that a single 
GeV photon pulse of 10 ys duration corresponds to a 
pulse energy of only about 100 pj, but to a power of 
10 TW. A single photon per pulse is in principle sufficient 
to reconstruct a non-trivial pulse shape [34j]. Decreasing 
the photon energy would enhance the photon yield, but 
would also increase the number of unwanted background 
photons. Increasing the collision energy would further 
enhance the number of photons produced, and could also 
increase the relative importance of the contribution of 
thermal photons compared to other kinds of photons [32l |. 

Concluding, we have studied the time evolution of the 
photon emission from the QGP. Since the QGP only ex- 
ists on the yoctosecond timescale, it naturally emits pho- 
tons at this timescale. We have shown that the emis- 
sion envelope can be influenced by the geometry, emis- 
sion angle, and internal dynamics like the isotropization 
time of the expanding QGP. For a particular parameter 
range, that is non-central collisions, large isotropization 
time, and an emission angle close to forward direction, a 
double-peak structure has been found within our model. 
Thus, pump-probe experiments at the GeV energy scale 
could be envisioned. New detection schemes would be 
required, and it should be explored whether tools and 
ideas from attosecond metrology 0, [35| can be scaled 
to zepto- or yoctosecond duration, utilizing physical pro- 
cesses in the GeV energy range, like electron-positron pair 
creation. Alternatively, determining the photon emission 
shape experimentally would give direct access to dynamic 
properties of the QGP, like its isotropization time. 

We thank K. Z. Hatsagortsyan and M. Strickland for 
helpful discussions. 
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